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ABSTRACT 

Context. The projection factor (p), which converts the radial velocity to pulsational velocity, is an important parameter in the Baade- 
Wesselink (BW) type analysis and distance scale work. The p-factor is either adopted as a constant or linearly depending on the 
logarithmic of pulsating periods. 

Aims. The aim of this work is to calibrate the p-factor if a Cepheid has both the BW distance and an independent distance measure- 
ment, and examine the p-factor for 6 Cephei - the prototype of classical Cepheids. 

Methods. We calibrated the p-factor for several Galactic Cepheids that have both the latest BW distances and independent distances 
either from Hipparcos parallaxes or main-sequence fitting distances to Cepheid-hosted stellar clusters. 

Results. Based on 25 Cepheids, the calibrated p-factor relation is consistent with latest p-factor relation in literature. The calibrated 
p-factor relation also indicates that this relation may not be linear and may exhibit an intrinsic scatter. We also examined the dis- 
crepancy of empirical p-factors for 6 Cephei, and found that the reasons for this discrepancy include the disagreement of angular 
diameters, the treatment of radial velocity data, and the phase interval adopted during the fitting procedure. Finally, we investigated 
the impact of the input p-factor in two BW methodologies for 6 Cephei, and found that different p-factors can be adopted in these 
BW methodologies and yet result in the same angular diameters. 

Key words, stars: variables: Cepheids — stars: individual: 6 Cephei — distance scale 
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1. Introduction 

The projection factor, or p-factor, converts the observed radial 
velocity to pulsational velocity and is a key parameter in the 
Baade-Wesselink (BW) type analysis and in distance scale appli- 
cations. Currently the p-factor in the literature is either adopted 
as a constant with the period (some of the values are given, for 
example, in Burki et al. 1982; Kovacs 2003; Groenewegen 2007 ; 
Feast et al. 2008 ) or involves a period-dependence in the form of 
p - a log(f) + b (the Pp relation, see, for example, Gieren et 
al. [19931 [2005: Nardetto et al. l2007ll20Tj9l Laney & Joner l2009l 
Storm et al. 1201 11 1. These p-factors are listed in TableQ] The lat- 
est derivation of a Pp relation by Storm et al. (2011. hereafter 
Sll), is based on the combination of two constraints: Galactic 
Cepheids possessing accurate parallaxes from the Hubble Space 
Telescope (HST, Benedict et al. 2002 2007) to constrain the in- 
tercept of the Pp relation, and Large Magellanic Cloud (LMC) 
Cepheids to determine the Pp slope by demanding that BW dis- 
tances to these LMC Cepheids be independent of the pulsation 
period. The p-factor is found to be insensitive to variations in 
chemical abundance (Nardetto et al. 1201 11 1. 

The Gaia mission will not present distances to Galactic 
Cepheids in the next few years (though it is scheduled to launch 
in 2013). Until then, the BW method is the best avenue for mea- 
suring the Cepheid parallax beyond HST and main-sequence fit- 
ting, but the method depends on the p-factor, which is known to 
depend on atmospheric physics, the circumstellar medium, and 
observational bias (e.g., the way radial velocities are measured). 



Table 1. Comparison of the p-factor. 



Reference 


a 


b 


Burki etal. (19821 






1.36 


Kovacs {20031 






1.35 


Groenewegen (120071 






1.27 ± 0.05 


Feast et al. 4200811 






1.23 ± 0.03 


Gieren etal. (|1993) 


-0 


03 


1.39 


Gieren et al. d2005t 


-0.15 


±0.02 


1.58 ± 0.02 


Nardetto et al. (2007 ( 


-0.075 


±0.031 


1.366 ±0.036 


Nardetto et al. (2009 ( 


-0.08 


±0.05 


1.31 ±0.06 


Laney & Joner ( 2009 ( 


-0.071 


± 0.020 


1.311 ±0.019 


Storm etal. d201ll 


-0.186 


±0.06 


1.550 ±0.04 



The goal of this paper is to extend the work of Sll. Because 
SI l's Pp relation was calibrated using Cepheids with HST par- 
allaxes (and LMC Cepheids), we demonstrate in Section 2 that 
the Pp relation can also be calibrated using other independent 
distances. The discrepancy of the p-factor for 5 Cephei is dis- 
cussed in Section 3, we explore in which way different methods 
affect the p-factor determination by using different angular di- 
ameter measurements and BW methods. The conclusion is given 
in Section 4. 
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Fig. 1. Top panel (a): Distance ratio as a function of BW dis- 
tances given in SI 1. Bottom panel (b): Distance ratio as a func- 
tion of pulsation periods in days. In both panels, the dashed lines 
indicate the case for And/ Aw = 1 and not the fit to the data. 
[See on-line edition for a color version.] 



Fig. 2. Calibrated p-factor for individual Cepheids using Eq. 
(1). The dashed line indicates the expected limit of the p-factor 
(p = 1.5, see text for more details). The Pp relation from Sll 
is included for comparison. Symbols are the same as in Fig. \T] 
[See on-line edition for a color version.] 



2. Calibration of the Pp relation using independent 
distances 

Because the p-factor is degenerate with the measured distance 
for a given Cepheid (for example, see Barnes et al. 2005a), the 
p-factor can be calibrated if a Cepheid has both the BW dis- 
tance (with an adopted p-factor) and an independent distance 
measured from other methods. That is 



Pnew = PBW X — , (1) 

"BW 

where D is the distance in parsec. This was performed, for ex- 
ample, in Groenewegen (2007 ), Feast et al. (120081 1. and Laney 
& Joner (2009), who compared the BW distances and geomet- 
ric parallaxes to calibrate the p-factor. Because HST parallaxes 
have been used to constraint the Pp relations in Sll, the de- 
rived BW distances are not fully independent of these parallaxes. 
Therefore, we adopted the revised Hipparcos parallaxes (van 
Leeuwen et al. 2007 ; their Tables 1 and 2) for nine Cepheids3, 
augmented with distances based on the main-sequence (MS) 
fitting technique from Turner (12010) for another 16 Cepheids, 
common to the Sll sampk^l to calibrate the p-factors. Note 
that the MS distance for TW Nor was updated based on the 
latest result from Majaess et al. (1201 1\ . The averaged distance 
ratio for these 25 Cepheids is < And/Aw >= 0.966, with 
a dispersion (cr) of 0.123. Figure Q] shows the distribution of 
And/ Aw as a function of BW distances (upper panel) and pul- 
sation periods (lower panel). A weak period dependence is found 

1 We excluded Y Sgr for the reason given in van Leeuwen et al. 
(2007 1, and RT Aur because it has a negative parallax. 

2 W Sgr was excluded from the sample based on the reasons given in 
Sll. 



for the distance ratio: And/Aw = -0.019(+0.065) logP + 
0.969(+0.067) with cr = 0.124, which is consistent with being 
period-independent. 

Distance ratios for individual Cepheids in the sample were 
used to calibrate the p-factors using Eq. (1). The resulting p- 
factors as a function of pulsation periods are presented in Fig. 
[2] where the error bars include propagated errors from both dis- 
tances and the Pp relation from Sll. A linear regression fit to 
all data yields p = -0.172(+0.086) logP + 1.462(+0.087), with 
cr = 0.107, and it is consistent with Sll's Pp relation. Figure 
[3] presents the correlation of the calibrated p-factor, p new , as a 
function of the distance ratios. Figures [2] and [3] reveal that about 
three Cepheids have calibrated p values (p ne w) that fall outside 
the expected limit of the p-factor: p > 1.5. These values indi- 
cate that limb brightening instead of limb darkening occurs in 
a Cepheid atmosphere (Sll). This can be seen from Eq. (6) in 
Nardetto et al. (2006a): p c - -0.18mv + 1.52, where p c is ge- 
ometric p-factor and uy is the limb darkening in V band. For a 
uniform limb darkening, uy = when p ~ 1.5, and hence p > 
1.5 implies a limb brightening (also, see Neilson et al. 120121 ). 
Since p new is degenerate with And /Aw, as shown from Eq. (1), 
Pnew > 1-5 for these Cepheids suggested that either the inde- 
pendent distances are overestimated or the B W distances are un- 
derestimated, or both distances are incorrect. Nevertheless, these 
Cepheids have p new within ~ cr p from the limit (where cr p is the 
estimated error on p ne w)- The discrepant point shown in Fig.|2]is 
FF Aql, with p = 1.89. This is due to the large difference in dis- 
tances from HST {jihst = 2.81 + 0.18 mas [milli-arcsecond], or 
D = 355.9+22.8 pc) and Hipparcos (7Thippai-cos = 2.05 +0.34 mas, 
or D = 487.8 + 80.9 pc). After removing FF Aql, the resulting 
Pp relation is: p = -0.159(+0.070) logP + 1.447(+0.070), with 
a dispersion of 0.064. The Pp relation from Sll would be a pre- 
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Fig. 3. Correlation of distance ratios and calibrated p-factors for 
Cepheids in the sample. The dashed line indicates the expected 
limit of p-factor (p = 1 .5, see text for more details). The slope of 
the data at the point where the intercept passes through the ori- 
gin, is 1.374, which is almost identical to -0.186 x < log(f) > 
+ 1.550 = 1.372, with a mean period of < log(P) >= 0.955. 
Symbols are the same as in Fig. Q] The outlier with largest error 
bars is FF Aql. [See on-line edition for a color version.] 

ferred relation, because it is calibrated with HST parallaxes and 
constrained from LMC Cepheids. 

Based on a slightly larger sample of Galactic Cepheids than 
that of SI 1, Fig. [2]reveals that the /^-factor relation may not be 
a constanfl Furthermore, the period dependency of Pp relation 
may not be linear, especially if the period-luminosity (P-L) re- 
lation is nonlinear (see, for example, Kanbur & Ngeow 2004, 
Sanda ge et al. 120041 Ngeow et al. [20051 |2009l Neilson et al. 
120101 and reference therein). From the perspective of the geo- 
metric p-factor, it depends on limb darkening, which in turn de- 
pends on the combination of luminosity, effective temperature, 
and gravity (Neilson & Lester 1201 11 . If the P-L and period- 
temperature (or color) relations are nonlinear, then the same 
holds for the geometric Pp relation. This nonlinearity could be 
in quadratic or in other formsQ. Nevertheless, confirmation or 
refutation of nonlinearity of the p-factor relation has to wait for 
accurate parallax measurements from Gaia mission for Cepheids 
in S 1 1 sample. Furthermore, based on spherically symmetric at- 
mosphere models, Neilson et al. (120121 1 found that the theoretical 
Pp relation is also nonlinear. 

Figure [2] also suggests a possible existence of intrinsic dis- 
persion on the Pp relation, albeit large error bars, which could 
be caused by systematics in the BW or independent distances, 
or both. Indeed, the /^-factor measured in this work includes all 
potential uncertainties in the BW method (or there is a prob- 

3 For completeness, the weighted mean for the p-factors shown in 
Fig.|2]is 1.306 ± 0.019 (<r = 0.187), or 1.304 ± 0.019 (cr = 0.148) if FF 
Aql is removed. 

4 For examples, two linear regressions with a break period at 1 days, 
or in the form of a constant + linear regressions, or a power law relation. 



Fig. 4. Comparison of the angular diameters for 6 Cephei as 
a function of pulsation phases from two methods. These two 
methods are the interferometric technique (Merand et al. 2005 
their Table 3) and IRSB technique using the prescription given 
in SI 1, with VK band photometric light curve data taken from 
Moffett & Barnes (TT9g4b and Barnes et al. (TT997l >. Epoch of 
2443674.144 days and period of 5.366316 days are adopted 
from Moffett & Barnes (1985). Error bars are omitted for clar- 
ity. The dashed curves are modeled angular diameters as a func- 
tion of pulsational phases based on the solution of Eq. (2), using 
D = 273.2 pc, see text for more details. 



lem in the implementation of the BW method). The dispersion 
may be a natural result of the width of the instability strip, i.e. 
period-dependent relations for Cepheids typically exhibit an in- 
trinsic dispersion (e.g., the period-color relation). Another pos- 
sible source for the dispersion may be dynamics in Cepheid at- 
mospheres. Assuming that pulsation amplitudes provide a mea- 
sure of atmospheric dynamics, then the significant dispersion in 
period-amplitude relations (as shown in Klagyivik & Szabados 
2009) suggests that dynamics can contribute to the significant 
dispersion of the Pp relation. Again, testing the existence of in- 
trinsic dispersion has to await parallax measurements from Gaia 
mission. 



3. The p-factor for 5 Cephei 

Merand et al. (2005) have directly measured the /^-factor, p - 
1.27 + 0.06, for 6 Cephei using interferometric measurements 
and the HST parallax. This value is consistent with the result de- 
rived from theoretical predictions (p = 1.27 + 0.01, Nardetto et 
al. 2004 ) or with the value of p — 1.24+0. 05 from the Pp relation 
provided in Nardetto et al. (|2009). The p-factor established by 
Merand et al. (2005) and Nardetto et al. (12009) is tied to radial 
velocities inferred from the cross-correlation method, whereas 
Nardetto et al. (2004) relied on applying Gaussian fits to deter- 
mine the radial velocities. In contrast, the p-factor derived from 
using the BW distance and HST parallax is 1.45 + 0.07 ( SI 1), 
about 2<x larger than the value given in Merand et al. (2005 1. The 



3 



Chow-Choong Ngeow et al.: Calibrating the projection factor for Galactic Cepheids (RN) 



predicted p-factor using the Pp relation from SI 1 is 1.41, which 
is still larger than the value from Merand et al. (2005 1. In this 
section, we examine the possible reasons for this discrepancy. 



3.1. Angular diameters ford Cephei 

Figure [4] compares the computed angular diameters as estab- 
lished from interferometric technique (Merand et al. 2005} and 
infrared surface brightness (IRSB) technique as a function of 
pulsation phases. It is clear from this figure that the angular di- 
ameters from these two techniques do not agree, especially at 
phases from ~ 0.2 to ~ 0.7, even though Kervella et al. (2004) 
have shown a good agreement for i Car. The mean angular diam- 
eters from interferometric and IRSB techniques are 1 .475 mas 
(with cr = 0.053) and 1.508 mas (with cr = 0.068), respectively. 
Cepheid angular diameters measured from interferometric ob- 
servations have been corrected using limb-darkening relations 
from model stellar atmospheres. Marengo et al. (12003) found 
that the model corrections introduced negligible errors and were 
largely insensitive to atmospheric dynamics. That analysis as- 
sumed plane-parallel radiative transfer, but Neilson & Lester 
(2011) found that for stars with low gravity (logg ~ 1-3), 
more realistic spherically symmetric model atmospheres pre- 
dicted greater variations of limb darkening than those found by 
Marengo et al. d20031 >. Neilson et al. (120121 1 also showed that 
angular diameter corrections from plane-parallel model atmo- 
sphere lead to an approximately 2% underestimate of the an- 
gular diameter, which gives a ~ 2% systematic underestimate 
of the p-factor. It is likely that the interferometric observations 
predict angular diameters that are smaller than predicted by the 
IRSB, which measured the angular diameter from the stellar 
flux and temperature-color relations. Nevertheless, Nardetto et 
al. (I2006bb showed that the derived distance to 5 Cephei is not 
affected by the limb-darkening variation with phases on the in- 
terferometric angular diameters. 

Another possibility is that the angular diameters from the 
IRSB method is overestimated. Neilson et al. (j2010) found that 
angular diameters measured using the IRSB techniques may be 
overestimated because of circumstellar media that cause an in- 
frared (IR) excess. This IR excess is presumably caused by mass- 
loss activity, evidence for which has been found recently on 6 
Cephei (Marengo et al. 120 101 Matthews et al. 1201 21) . In particu- 
lar, the /T-band flux excess for 5 Cephei is 1.5% (Merand et al. 
2006), which translates into a difference of AK = 0.016 magni- 
tudes. Using the surface brightness relation (Sll), this flux ex- 
cess translates into a difference in angular diameter, A log 6 = 
0.004. For a mean angular diameter of ~ 1 .48 mas, the IRSB 
technique would overestimate the angular diameter by AG = 0.01 
- 0.02 mas or AO/0 — 1%. /if -band flux excess may explain a 
significant fraction of the difference seen in Fig. [4] However, a 
1% difference in angular diameter leads to a 1% difference in 
p-factor, suggesting the A'-band flux excess is insufficient to ex- 
plain the difference in p-factor for 6 Cephei found by Merand et 
al. (120051 and Sll. 



3.2. The p-factor derived from the BW method 

In the BW analysis, the angular diameters as shown in Fig.|4]can 
be modeled according to the following equation: 



as 



•Barnes et al. (2005) 
oStorm et al. (2004) 
"Bersier et al. (1994) 




phase 

Fig. 5. Folded radial velocities for 6 Cephei, using an epoch 
of 2443674.144 days and period of 5.366316 days. The data 
were fitted with an 8 th order Fourier expansion, shown by the 
solid curve. Data from Bersier et al. (1994) needed to be shifted 
slightly to minimize the scatter in the combined radial velocity 
curve. Error bars are omitted for clarity. [See on-line edition for 
a color version.] 

Table 2. Derived p-factors for S Cep. 



Data source for 6 


Include all phases 


Exclude phases > 0.8 


Merand et al. (120051 
Storm etal. (12011* 


1.32 ±0.05 
1.50 + 0.05 


1.26 ± 0.08 
1.58 ±0.08 



IPC f 

= s -p— J o [v 8 («-r]* 



(2) 



where is the pulsational phase, Vr is the radial velocity curve, y 
is the systemic velocity, and C - 0.57749 is the conversion fac- 
tor (for P in days, D in parsec, velocities in km/s and 9 in mas). 
Given an adopted p-factor, Eq. (2) can be used to derive the dis- 
tance (and radius) to a Cepheid. Conversely, the p-factor can be 
determined if the distance is known a priori. The radial veloc- 
ity data for 5 Cephei are taken from Bersier et al. (11994) . Storm 
et al. ( 120041) . and Barnes et al. (2005B. The combined radial 
velocity curve is fitted with an 8 th order Fourier expansion, as 
presented in Fig. [5] and a distance of D — 273.2 pc was adopted 
for 6 Cephei (Benedict et al. 2002). The fitting procedure was 
performed using an OSL bi-sector algorithm from SLOPE (Isobe 
et al l 19901 ). The fitted p-factors, based on Merand et al. (120051) 
and Sll angular diameters, are summarized in Table|2] The solu- 
tions (9q and p) were subsequently used to construct the angular 
diameter curves using Eq. (2). These curves are shown in Fig. [4] 
and agree well with the empirical angular diameters. 

The derived p-factors using all data points, p = 1.32 + 0.05 
and 1 .50 + 0.05, are consistent with p = 1 .27 from Merand et al. 
(120051 ) and p — 1 .45 from Sll, respectively, although the derived 
values are higher by 0.05. The p-factors from using the Merand 
et al. (120051 ) and Sll angular diameter data disagree at ~ 2.5<x 
level, due to the disagreement of the angular diameters shown in 
Fig. |4] The difference of the p-factor stems from the manner in 
which the radial velocities were treated, since a similar (or the 



4 



Chow-Choong Ngeow et al.: Calibrating the projection factor for Galactic Cepheids (RN) 



same) period and distance for 5 Cephei were adopted in Merand 
et al. (2005), S 1 1, and this work. Merand et al . (120051 considered 
radial velocity data either from Bersier et al. (11994) or Barnes et 
al. d2005bl l: while Sll included several additional data sources 
to those adopted in this work. For fitting the radial velocity data, 
Merand et al. (2005 ) applied a four-knot periodic cubic spline 
interpolation, in contrast to S 1 1 and this work, where a Fourier 
expansion was used to fit the radial velocity datcQ 

Finally, Merand et al. (120051 1 used all data points during the 
fitting procedure, while Sll excluded data points with <p > 0.8 
due to the deviation of angular diameters at these phases during 
the fitting (probably caused by shock waves at minimum radius, 
see also Fouque et al. [2003, Kervella et al. I2004I I. In Table |2l 
we also include the fitted p-factors if data points for <p > 0.8 
are excluded, as in S 1 1 . These p-factors agree with the p-factors 
derived from using all data, but they disagree at ~ 2.8<x level 
(the reason is again owing to the disagreement of the angular 
diameters). Nevertheless, it is clear from Table|2]that the adopted 
phase intervals will affect the fitted p-factor. 

3.3. p -factor from other BW methodologies 



Feast et al. ([2008) and Laney & Joner (2009) employed a dif- 
ferent BW method to derive the distance to 6 Cephei. Their 
method is described in more detail in Laney & Stobie (1995 
hereafter LS95). After deriving the BW distance, the p-factor is 
recalibrated using geometrical parallaxes. We focus on the re- 
calibrated p-factor from Laney & Joner (2009l, because they 
adopted the same HST parallax as in Sll, which yields 1.289 + 
0.06 1 . Their p-factor is also ~ 2cr smaller than the p-factor from 
Sll. On the oth er hand , the BW distance of 268.8 + 7 pc from 
Laney & Joner (120091 1 is in excellent agreement with the BW 
distance of 266.7 + 5 pc from Sll, even though the adopted p- 
factor is different: 1.27 versus 1.41. This reflects that the differ- 
ence in the adopted p-factor may compensate for the difference 
in assumptions and methodologies in these two BW methods. 

This can also be seen from the following. In Sll, the pho- 
tometric surface brightness relation can be expressed as Fy = 
a(V - K) + p and F v = 4.2207 - 0.1 V Q - 0.5 log 0((f>). Equating 
these two expressions and re-arranging them, yields 



V = a'(V-K) -5\o g e(cf>)+p', 



(3) 



where a' = -10a and ft' = 10(4.2207 - 0) are constants, with 
a = -0.1336 and f3 = 3.9530 as given in Sll. In the LS95 
methodology, the radius in V = a(V - K) Q - 5 log(R + AR) + b 
can be converted into angular diameters: 



2.2 
1.6 



•Fitting a, b' 
D 0nly fitting 



and 
S„ in 



8 in equation 
equation (3) 



(4) 



□ □□□□□□□□□□□□□□□□ D- g DDnnnDDDnnn 



1.3 1.3 1.4 1.5 

Input p-factor 

Fig. 6. Fitted values of b' (upper panel) and 9o (lower panel) in 
Eq. (4) as a function of input p-factors. Filled circles are the 
fitted results using Eq. (4), and open squares are the fitted 6q 
using Eq. (3). Note that at 273.2 pc, various published period- 
radius relations give 6q in the range from ~ 1 .4 mas to ~ 1 .7 mas 
for 6 Cephei. 



of input p-factors (in contrast to 1.336 given in Eq. [3]). The 
fitted coefficients give the same sum of residuals squared for all 
input p-factors. For comparison, we also fit the 6>o in Eq. (3) for 
the same input p-factors. The fitted results are shown as open 
squares in the lower panel of Fig. [6] This shows that in the LS95 
methodology, 0q is more sensitive to the input p-factor than the 
Sll methodology. Note that if the mean radius of 6 Cephei is 
41.9.Ro (Feast et al. 1200811, then O = 1.45 mas at the distance 
found in Laney & Joner (2009] D — 268.8 pc). Based on the 
lower panel of Fig. [6] the required input p-factor to obtain 6*o = 
1.45 mas is 1.38 and 1.48 for the LS95 and Sll methodology, 
respectively. Even though the p-factors found here are higher 
than 1.27 and 1.41, this example demonstrates that different p- 
factors adopted in the LS95 and Sll BW methodology result in 
the same angular diameter. Hence, it is possible to obtain similar 
distances with different input p-factors in different BW methods 
by compensating for the fitted coefficients. 



V = a(V -K) -5 log 0(cf>) + b', (4) 

where b' = b — 5 log(Z)/2) (here, D is treated as a nuisance pa- 
rameter). In both Eq. (3) & (4), the 9{<p) is given in Eq. (2). In 
contrast to the Sll methodology, the coefficient a and b' , as well 
as Oq, in Eq. (4) are solved (or fitted) and not fixed as constants. 

Using the same VK photometric data as in Fig. [4] and the 
radial velocity curves as presented in Fig. [5] we fit coefficients 
a, b' and 6>o in Eq. (4) for a range of input p-factors (from 1 .20 
to 1.50 with a step size of 0.01). The results are presented in 
Fig. [6] for V and 0q, the fitted value for a is 1.292 regardless 

5 The difference of p-factors from Sll and Table |2]may include the 
different order of the Fourier expansion: Sll adopted a lower order fit 
while in this work higher order terms are included. 



4. Conclusion 

We have calibrated the p-factors for Galactic Cepheids using the 
simple fact that the p-factor is degenerate with distance, by using 
Cepheids with both BW distance and an independent distance 
measurement. The latest BW distances for Galactic Cepheids 
are provided in Sll, and the independent distances were ei- 
ther adopted from Hipparcos parallax measurements or from the 
MS fitting. We also compared the p-factors established for 6 
Cephei by Merand et al. (2005 interferometric technique) and 
Sll (IRSB). Several physical aspects might explain why the 
aforementioned groups disagree concerning the p-factors: (a) in- 
consistent angular diameters as inferred from the interferometric 
and IRSB techniques; (b) a difference in the adopted radial ve- 
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locity data and the interpolating/fitting of these data; and (c) a 
difference in the sample selected to fit the data (i.e. excluding 
data with (f> > 0.8). We also examined a different BW method- 
ology than SI 1, as presented in LS95 and Feast et al. d2008), 
which adopted a different p-factor and yet obtained almost the 
same distance to 5 Cephei as in S 1 1 . This is because the surface 
brightness coefficients are fitted from data in LS95 methodology, 
which can compensate for a different input p-factor. 

Even though there are currently only ten Cepheids with ac- 
curate geometrical distances measured from HST, a larger num- 
ber of Cepheids will have better parallax measurements after the 
launch of Gaia. These parallaxes can be used to verify the Pp 
relation derived in SI 1, and to examine the nonlinearity and 
possible intrinsic dispersion of the relation. Furthermore, di- 
rect p-factor measurements are being obtained for a few addi- 
tional Cepheids using interferometric techniques (Kervella 201 1 
— private communication), and the results will be employed to 
evaluate the viability of the conclusions presented here. 
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